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I. INTRODUCTION 

Recently, the so called ridge phenomenon has been ob- 
served at RHIC [l[ (particle yield enhancement in narrow 
azimuthal angle window A<f) but broad pseudo-rapidity 
window A 77). It is a jet- related long range correlation in 
rapidity. For us it is particularly interesting that the so- 
called soft ridge phenomenon, a similar longitudinal ra- 
pidity correlation but without jet trigger, is found in col- 
lisions of higher centrality events [1, H[ . Although there 
is no theoretical model that can quantitatively reproduce 
all experimental data related to this phenomena, several 
models have been proposed aiming at the explanation of 
ridge data 0-[lo|. 

Actually, correlations over several rapidity units can 
only occur at the earliest stages of heavy-ion colli- 
sion, when different rapidity regions are still causally 
connected [H, @. It means that the origin of this 
phenomenon must be placed almost instantaneously 
after the collision of two nuclei 0, [|| , namely before or 
around the initial stage of the hydrodynamical evolution 
of the fireball created by the relativistic heavy-ion 
collisions. On the other hand, observation of correlations 
is related to period after freeze out, what means that 
the long-range correlations can survive the whole stage 
of hydrodynamical evolution of the matter. Since little 
attention has been given to this problem so far it is 
therefore naturally to ask what are the influence of the 
long-range correlations on the fluid dynamical evolution 
because, this will be the subject we shall investigate 
here. Because such dynamical long-range correlations 
may cause weak dissipation and viscosity in the fluid 
we propose a kind of extreme model for the relativistic 
dissipative hydrodynamics with small dissipative and 
viscous terms, which originate entirely from the long- 
range correlations among the fluid elements. 

A few words on hydrodynamical models used in 
particle production processes are in order here. The 
ideal fluid dynamics [llT - [l5j have successfully reproduced 
experimental results obtained at Relativistic Heavy-Ion 
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collider (RHIC) (for a review, see for example, fl7l|). 
This success brought us a new understanding of the 
matter created at the RHIC. It turned out that this 
matter behaves more like an ideal fluid rather than the 
weakly interacting parton plasma fl8j . It is thought 
that this is because the mean free path of the particles 
composing it is small comparing to the system size or 
to the typical scale in the fluid dynami cs ( strongly in- 
teracting quark gluon plasma, sQGP |19l - [2ll |). However, 
there are noticeably differences between predictions 
from the ideal fluid model and the cor resp onding data 
for the flow parameters obtained in [22l424j [Hj]. Hence, 
some revisions to the ideal fluid picture turned out to 
be necessary and dissipation and viscosity effects have 
been introduced to the perfect fluid. This caused some 
problems with the proper formulation of the relativistic 
dissipative fluid model (25rl33j . In a simple relativistic 
extension of the Naiver-Stokes equation (so-called first 
order theories 03 HH), it is known that the causality 
is violated [3rj l37l| and the solution of the equation is 
not stable against small perturbations [H, 39] . Over the 
past few years a considerable number of studies have 
been made aiming at the solution of these problems (see, 
for example Ref . g3"Sl ) • 

The purpose of this paper is to formulate the relativis- 
tic imperfect hydrodynamics including the effect of long- 
range correlations and to find its connections with the 
usual dissipative hydrodynamics. Note that, in the pres- 
ence of the long-range correlations over the thermal equi- 
librium scale, the microscopic state of the particle that 
composes the fluid does not obey the usual Boltzmann- 
Gibbs statistics any longer 43]. Also thermodynamic 
functions such as the energy density and pressure among 
the separate fluid cells will be affected. The state in the 
presence of the long-range correlations, which is gener- 
ally different from the usual equilibrium state (it is a 
kind of a non-equilibrium state), can be the so-called 
stationary state near the equilibrium. It approaches to 
the equilibrium state with a finite (or infinitely long) re- 
laxation time under a non-equilibrium thermodynamical 
constraint which may make a kind of narrow valley of 
the free energy generated by dynamical correlations [44] . 
This relaxation process depends on the initial state se- 
lected in the valley of the free energy. Such dependence 
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on the initial state for the relaxation processes can be 
then regarded as a kind of memory effect induced by the 
long-range correlations. 

The relaxation process of the stationary state is 
called pre-thermalization [44| (see also [45j]) and leads 
to Tsallis's statistics [If|. It is usually considered to be 
important for describing correlated system and appears 
when the phase space is reduced by the correlations 
among constituents of the system. In [47| it is shown 
that there is a certain correspondence between the fluid 
dynamics based on the Tsallis statistics and the dissipa- 
tive fluid dynamics based on the usual Boltzmann-Gibbs 
statistics. Thus, the dynamical long-range correlations 
may result in some dissipative contribution in the 
relativistic fluid. 

It should be stressed at this point that situation 
considered here is rather special. Our model assumes 
that, because of action of the long-range correlations, 
some stationary, near equilibrium, state of fluid is 
formed replacing the usual equilibrium state. The 
natural question is how real is such assumption. So far 
there is no hard evidence for creation of such state. On 
the other hand, there is growing evidence of a power 
law-like behavior of the transverse momentum spectra 
[13, HI] in the mid p T region (2-6 GeV/c) and this 
may indicate the existence of such stationary state 
[47j . The necessary condition for the existence of such 
pre-thermalization state in the matter created in the 
relativistic heavy-ion collisions (for example at RHIC) is 
that its relaxation (or equilibration) time is long enough 
comparing with the typical hydrodynamical time scale, 
such as L hydro /c s , where L hydro is some typical length 
scale of fluid and c s is sound velocity. All these calls 
for some very detail considerations which are, however, 
outside of the limited scope of the present paper. Here 
we consider therefore only a kind of thought experiment 
in which appearance of the stationary state due to the 
long range correlations is assumed and its consequences 
investigated in detail. 

The long-range correlation influence not only the 
microscopic state but also the macroscopic state of 
the fluid such as velocity vectors. Therefore, instead 
of specifying the type of the long-range correlations, 
we shall assume that all information on them are 
included in a flow velocity vector field, u M (x), and in 
the non-equilibrium thermodynamic functions like local 
energy density, e(x), and pressure, P(x) (here x denotes 
a space-time four vector). 

To quantitatively investigate the effects of the long- 
range correlations (LRC) on the relativistic fluid, we 
shall propose a thought experiment assuming that it is 
possible to switch on and off the LRC for a perfect fluid 
in such a way as not to break the causality. If one can 
turn on the switch of the LRC that acts on the perfect 
fluid at a proper time x° = r on , then, after r on , the 



fluid starts to expand as a dissipative fluid. We assume 
that solution of the perfect fluid, which is obtained 
in the case without LRC, is known for all space-time 
points (i.e., by solving ideal fluid dynamics after r = r on 
we know the energy density, e cq (x), conserved charge 
density, n cq (x), flow velocity, U ll (x), and equation of 
state P eq (x) — P eq (e eq , n eq ) for all space-time point x 
[49l IHoj . It is then assumed that to know evolution of 
changes (defined as A(x) = e — e oq , LT(x) = P — P cq , 
8n(x) = n — n cq and <5u M = u M (x) — U fJ '(x)) is equivalent 
to solving equations of resulting dissipative hydrody- 
namics. As will be seen below, since A, LI and Sn can be 
expressed as a function of the enthalpy, h cq = e cq + P cq , 
n oq , and Su^, in order to investigate the dissipative 
effects of the LRC is sufficient to know only the evolution 
of Su^^x). Hereafter, we shall call the set of solutions 
of the ideal hydrodynamical equations, e cq (x), P oq (x), 
n cq (x) and U y '{x), the reference fluid fields (RFF) for 
the corresponding dissipative fluid field. One may also 
define RFF by switching-off the LRC for a given fluid 
with dissipation that originated from LRC. 

The paper is organized as follows. In Sec. II, we briefly 
review the relativistic fluid dynamics. Then, in Sec. Ill, 
we introduce the RFF which satisfy ideal hydrodynamics 
and we express the energy-momentum tensor in terms 
of Su^ and RFF. In Sec. IV equation on <5it M is found 
which agrees with the relativistic dissipative hydrody- 
namic equations. Sec. V contains the off-equilibrium 
entropy current corresponding to our model. We close 
with Sec. VI containing the summary and some further 
discussion. 



II. THE RELATIVISTIC FLUID DYNAMICS 

Basic equations of the relativistic hydrodynamics for 
arbitrary dissipative fluid consist of by local conservation 
laws for the energy-momentum and conserved charges. In 
addition to these local conservation laws, the second law 
of thermodynamics is required. These equations are gen- 
erally expressed by the covariant derivatives (55| of the 
energy-momentum tensor, T llv {x), the conserved charge 
current, N IJ, (x), and the entropy current S^(x); 

= 0, (la) 
tf& = 0, (lb) 
> 0. (lc) 

The general form of the energy-momentum tensor used 
here is given by 

= eu^u v - PA^(u) + 2W^u v) + tt"" , (2) 

where e = T^u^Uy is energy density, P = 
-iA^w)^ is pressure, W* = u l/ T uX A^(u) is energy 
flow and tt^ = T^") is shear stress tensor. The fluid is 
moving with four- velocity u^lx) such that u^u^ = 1. We 
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introduce transverse projection tensor A '"'(it) = — 
u^u v with general metric for the space-time consid- 
ered. This tensor has the property that u^A^ v (u) = 
u u A^{u) = 0. The notation used is: A&B") = 
^(A^B" + A v ' B^) denotes symmetric tensor defined by 
two four-vectors A^ and B v and the angular bracket no- 



tation A^h 



A a ? 



rep- 



resents symmetric and trace free part of tensor A 11 " . 

The net charge and entropy currents are generally 
given by, respectively, 



S" 



(3) 
(4) 



where $ M = A^(u)S v is entropy flux. In this paper we 
shall use the so-called Eckart frame [3J] and fix arbitrary 



of choice of frame for the four flow vector 



In what 



follows we assume that baryon number, iV, is conserved. 

If the local thermal equilibrium is achieved (or, if the 
LRC is switched off in our thought experiment and af- 
ter the relaxation process following it has been com- 
pleted), the entropy flux <I>' 1 should also disappear, its 
density reaches the maximum, and it is locally conserved. 
The energy flow vector W M and viscous shear tensor it'- 11 ' 
should also disappear. One gets therefore a perfect fluid, 
corresponding to Eq.@, © and ((U), for which one has 



N? q = n cq (I»[/'\ 
= Scq (T, M )C/^, 



Peq(£e q ,n cq )A^([/), 



(5a) 
(5b) 
(5c) 



where T = T{x) and fi = fJ-(x) are the tempera- 
ture and (baryonic) chemical potential field given by 
£ e q(x) and n e q(x). The equation of state P cq (x) — 
Peq{£eq(x),n cq (x)) for the perfect fluid are assumed to 
be known. Here, U^ix) is the flow vector appearing af- 
ter relaxation has been completed. 



III. RELATIVISTIC HYDRODYNAMICS IN 
THE PRESENCE OF LONG-RANGE 
CORRELATIONS 

A. Reference fluid field (RFF) 



charge density, n C q(x), pressure P cq (x) and three un- 
known components in the fluid velocity, U^ix). On the 
other hand, there are 6 equations: the first two equations 
of Eq.© and equation of state, P cq (x) — P cq (e,n cq ) 
(equation Eq. (|6cl) . the local entropy conservations, can 
be derived from the previous two equations.). 



B. Effects caused by the LRC 

After switching on the LRC in our thought experiment 
(i.e., for a non-equilibrium dissipative situation) we have 
to solve Eq.(U) with Eqs.©, © and (@J. The pres- 
sure P is different from the isotropic equilibrium pres- 
sure, P eq (e eq , n cq ), which can be obtained by equation of 
state (EoS) via the equilibrium energy density, e cq , and 
the baryon number density, n oq (i.e., P ^ P e q(e e q,n e q)). 
The difference P — P eq is usually regarded as a bulk pres- 
sure, II = P — P eq . For energy density it is also natural 
to introduce bulk energy density, A = e — e cq , because 
the microscopic occupation probability does not follow 
the Boltzmann-Gibbs statistics any longer. Hence, the 
energy density e in Eq.([2]) cannot be expressed as a func- 
tion of the usual temperature T and (baryonic) chemical 
potential /i such like £ C q(T, fx). Therefore, one can write 
thermodynamical quantities in the presence of the long- 
range correlations, 



£ 

n 



- cq 



(7» + A, 
n eq (T, (j.) + 8n 



P -Pq(£eq? ^eq) ~\~ H. 



(7a) 
(7b) 
(7c) 



Now, the energy-momentum tensor Eq.Q can be de- 
composed by the flow vector U fi (x) introduced at the end 
of Sec. II, which is generally different from u M (x). De- 
noting by 5u^ the difference between u^(x) and U^(x), 
one can write 



u^{x) = U^{x) +5u fl (x). 



(8) 



Since both u' 1 and C/' 1 must satisfy normalization con- 
ditions as physical flow vectors, i.e., u^u^ = I and 
U^Uu = I, they are related in the following way: 



Let us introduce reference fluid field (RFF), by which 
we understand a set of solutions for a given perfect fluid 
without the LRC (because LRC is assumed to be the 
origin of dissipation and viscosity, one has perfect fluid 
if LRC are not present). In our case the perfect fluid 
dynamics is given by specifying 



rpflf 

1 cq-fj, — 

N>* = 

cq;/j, 

CM - - 



o. 
o. 

0. 



(6a) 
(6b) 
(6c) 



Note that there are 6 unknown variables in the perfect 
fluid hydrodynamics: energy density e cq (a;), conserved 



2 M 



u^5u a = 



7 



(9) 



It is desirable to describe the number of unknowns 
before moving to the main task. In dissipative hydro- 
dynamical model there are 14 unknowns: A, n, Sn, 3 
components of (because condition W^u^ — 0), 5 
components tt^ (symmetric and traceless and satisfying 
conditions ir^u^ — 0) and 3 components of Su^. On the 
other hand, e cq , n cq and 3 components of the [/ M (i.e., 
5 variables) and EoS, P eq = P e q(e e q,n eq ), are known for 
the reference fields. As will be shown in Sec. Ill D and 
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IV C, the 3 unknowns, A, II, Sn, and the 8 unknowns in 
and 7r Ml/ , total 11 unknowns are function of the 3 un- 
known of 8u^ and 6 known variables describing the RFF. 
Note that the number of unknown three variables for <5w M 
coincides with the number of spatial dimensions. They 
depend on the nature of the LRC (for example, if the cor- 
relations are totally spatially isotropic for a rest frame 
of U 11 , the number of unknowns is reduced to unity.). 
Therefore, we need to derive one equation from Eq. (|lap 
and Eq.flTET) (besides eqs. foaj) , (|6bf and EoS) to deter- 
mine one unknown. This is, as shown later, Eq. (|30ap or 
equivalently Eq. (|3"Tj) . 



C. Tensor decomposition by means of the RFF 

Introducing projection operator accompanying U^, 
A» v (U){= A"") = - UW, one can decompose 
Eq.© as 

T »u = i jjujjv _ p A ^([/) + 2W^ l U v) + tt^, (10) 
where e and P are 



rpafi 



T ap U a Up=e + 3Z, 



P 



1 



T af) A al3 (U) = P + E, 



(11a) 
(lib) 



respectively. The quantity E that appears both in the 
Eq. (|lla[) and (jllb|) is explicitly given by 



E = h(j 2 - 2 7 ) - - 7 ) + p' 1 



(12) 



where h = e + P is natural extension of the equilibrium 
enthalpy h cq = e cq + P oq and 



i = WJu", 7" = n ll Ju f *Su 1 '. 



(13) 



For the corresponding energy flow and shear viscosity 
tensor, we have 

= (1 - l){hX^ + Y^) - - iX» - 7"^, (14a) 

|(A^A^ + A^)-iA^A 



-liv 



+ hX»X v + 2X { ^Y v) + 2U^ l Z u \ 

respectively, where 

X» = A£(U)u a =Su^ + jU^, 
= A£{U)W a = W+j'U", 

= TT^SUu. 



(14b) 



(15a) 
(15b) 
(15c) 



which are vectors that can be constructed from <5u M , 
and ir^ . One can write then that 



Here we approach the essence of our model. In the case 
when e = s cq (x) and P = P cq (x) one has that 



A 

n 



-3E, 
-E. 



Assuming now that [5f 



= 0, 



(17a) 
(17b) 



(17c) 
(17d) 



one can link T^ v to its equilibrium correspondence 
T£ q . Physically, when the LRC is switched off the flow 
velocity (and also the related thermodynamical quan- 
tities) can not change into corresponding of the RFF 
instantaneously. The flow vector field u^{x) approaches 
its equilibrium limit U li {x) in a finite proper time. One 
may regard such change as a kind of relaxation process. 
In the next Section ITVl we shall derive the corresponding 
evolution of Su^ . 

Conditions Eqs. pT|) can be also mathematically ex- 
pressed by the following transformation: 



(18) 



Assumption ()18[) is natural because dissipative fluid 
must correspond to some ideal fluid, which can be 
obtained by switching off the LRC. It is our basic 
assumption deciding on the form of the energy flow W 1 * 
and shear viscosity tt^ , as will be shown below. 



D. The model energy-momentum tensor 

The assumption discussed in the previous Section con- 
cerning the energy-momentum tensor T^ v determines 
also the complete form not only of II but also both of 
W and tt^. From Eqs.(T7a} and ([i~7bl) . we obtain that 

M 



i.e., that 



A = 311, 



e = £ cq (T» + 3II, 

p = p eq (T, M ) + n. 



(19) 



(20a) 
(20b) 



A + 3E, P = P e < 



n 



(16a) 



Proceeding further, conditions = and = 
determine the explicit form of the energy flow vector 
and the shear tensor in the energy-momentum tensor 
in Eq.(fT]). To give the W 1 and tt^ correct physical 
meaning of the, respectively, energy flow vector and 
shear viscous tensor, we require that W^u^ = (or 
W^Suf, = 7') and that tt^u v = (or -K^du^ = Z u ). 
The later condition can fix the form of the four vector 
Z>> = - [(1 - 7)E + 7'] X» - i'U^. The former condi- 
tion W^u^ = gives a relation between E and 7'. From 
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the definition 7" = Z^Su^, eliminating and 7', we 
can also find relation between 7" and E. When these 
expressions for 7' and 7" are put into Eq.(fl2|). then one 
finds that 



7(2 - 7) 



(2 7 -l)(2 7 -3) 



(21) 



Eliminating 7', 7" and Z^ 1 from the equation = 0, 
one can obtain expression for in terms of U^ 1 and 
8u^. This result for put into the definition of 
leads to r^ 1 = -{h + Y)X». Substitution Y*, Z^ 
and 7" into the equation tt^ = 0, one gets expression for 
7r Miy . Finally, one finds that the energy flow vector and 
the shear viscous tensor have the following form: 



7T 



- [jh 4 - 3E] C/ M , 



(22a) 



= 7 2 /i 4 CW -EA^([/) 
+ h 2 Su t "Su u - 2 [2S - 7/14] tfu^E/") , (22b) 

where h n —i 2 4 = A + nE. Since h = s + P, one obtains 
that h = h cq + 411 (see Eq. ([20a]) and Eq. (l20bj) ). where 
/i cq = e cq + P cq . Using this relation and Eq. (|2"Tl) . one 
can write n(= — E) as function of 7 and of equilibrium 
enthalpy h cq as [581 ] 



n = -i 7 (2-7)ft eq (r,M)- 



(23) 



By using Eq. (j2"3")l . one can rewrite Eqs. (|22al) and (I22b[) in 

more compact form, such as 

W" = fc eq ( 1 -7)v M . (24a) 
tt^ = h^tpf +nA^"(u) = h cq 5u^8u v \ (24b) 

where 

^ = 7^ - (1 - 7 )<W. (25) 

Hence, our model energy-momentum tensor is finally 
given by 

= [e cq + 3U}u^u u -[P cq + U]A^(u) 

+2h cq [l - j](p^U^ + heqSuHu^ . (26) 

Since reference fields (e oq , n eq , P eq and U 11 ) are assumed 
to be known, the number of unknowns in this model 
energy-momentum tensor, Eq. (|26[) . is only three (let us 
stress here that this is because of the Eq. (fT8]) ). 

It should be noted that expressions Eqs. (|24a|) and 
(|24bl) . one finds tensor relations 



-3n (1 - 7 ) 2 /i eq , 



•k^Wu = -2IW, 



6ir, 



(27a) 
(27b) 
(27c) 



which are same as found in the Non-extensive/dissipative 
correspondence [47| providing 

&eq[l - if W q [l + 7 9 ] 2 , 



IV. EVOLUTION OF Su"(x) 

From the above discussion it is obvious that the most 
important quantity to be discussed in detail is 8u^{x). 
We shall therefore derive here equation which (5it M should 
satisfy. 



A. Dissipation caused by the LRC 

Let us start with effects of dissipation caused by 
switching on the LRC. Hereafter, we use , for simplicity, 
metric with gf 1 " = 0. The relativistic hydrodynamical 
equation Eq. ljlaj) can be always decomposed by contrac- 
tion with it M and A x (u) in the following way: 

u^e +{e + P)u£ 

+^ + M«"^ + 7^r]=0, (28a) 
[(e + PK + - d^PA x ^{u) 

+W x u^ + [u»W.l + n^]A x (u) = 0. (28b) 

One can also write the corresponding hydrodynamical 
equations for the RFF: 

C/^e cq + (e cq + P eq )C/£ = 0, (29a) 
(e cq + P cq )U»U^ - A x »(U)d^P cq = 0. (29b) 



Inserting Eqs.flHJ into Eqs. (|28al) and (I28b|) and then sub- 
tracting from them, respectively, Eq. (|29ap and (|29bl) one 
obtains 



lK q;fl + d^M 1 + * - 0, 



(30a) 



K^V X + d„P eq [2U^5u x) + Su^Su x ] + ^u x = 0, (30b) 

where h% q = h cq U^ is the enthalpy current in the equi- 
librium state and = h^UuSuY^. By eliminating ^ from 
Eq. (|30b|) using Eq. pOaj) one finds that Eq. (|30b)l is equiv- 
alent to Eq. (|29al) . This means that Eq. (|30b[) is not in- 
dependent equation. The unique independent equation 
obtained from Tfff — with model energy-momentum 
tensor (f2"6")l (besides Eqs. (|29a[) and (J29b|) remains there- 
fore Eq. (l30ap . which can be expressed as a constitutive- 
like equation: 

n + rnri = ^(l-7){ 7 e cq + V A1 P cq( S U ' 1 



+/i eq ( 7 6 + C/^)} r n , (31 



where X = U^X.^, VX = A" V (U)X. V and 9 = Ufa. 
The rn, which resembles a relaxation time, is given by 



-h eq (x)/h cq (x) 



(32) 



where 7 9 denotes the 7 that appear in 47 1 



and it is determined by the property of the reference 
field h e q(x), i.e., by the equilibrium enthalpy. Because 
Id is function of <5it M , the Eq. (|3"Tj) presents therefore the 
condition that <5u M should satisfy (as was also the case 
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for Eq.(3Dg}). 

Let us estimate rn roughly, starting, for simplicity, 
from 1+1 dimensional Bjorken scaling solution. The en- 
thalpy evolution is then given by the /i eq (r) ~ (to/t) 4 / 3 , 
where r is some proper time and tq is initial proper time 
needed to equilibrate the fluid without any long-range 
correlations. In this simple case, one finds that Eq. (j3"!2)) 
results in rn = §r. It means that, if the fluid at the initial 
stage is slightly different from the local equilibrium, it is 
difficult to reach its thermalization as proper time passes. 

In the more realistic 1+3 dimensional case, since the 
time derivative of the enthalpy \h cq \ in Eq. (|32[) must be 
much lager than this simple estimation, the rn should be 
smaller than the previous estimation. However, in this 
case the freeze out time, Tf D , is also shorter than in the 
1+1 dimensional case. Therefore, whether the stationary 
state (which is assumed to be create at the initial stage 
in the relativistic heavy-ion collisions) can arrive to its 
equilibrium depends on both the relaxation time and the 
freeze-out time of the created matter. Note also that the 
relaxation time rn does not depend on 5u^. The rn, in 
our model, is determined only by the nature of the equi- 
librium state. On the other hand, the equilibrium time 
r eq should depend on the Su^, which shows how much 
away we are from the equilibrium, and it is obtained by 
solving Eq.(l3Tj). 



B. A possible model for LRC 

We shall now discuss a possible model for the LRC 
in the hydrodynamic approach. Notice that we have 
only three unknown variables: Su , 5u 2 and 5u 3 . This 
is because A, LI, and tt^ are expressed in terms of 
<5w M and RFF whereas 6u° can be determined by the 
Eq.©. In Section IIIC it is shown that, because of the 
charge conservation, Sn can be also written in terms 
of 8u^ . Assuming that the LRC are spatially isotropic 
5u^ should also be isotropic in the rest frame of RFF 
velocity 

Let us consider a space-time point and the local rest 
frame for the reference flow vector XJ^ix) — (1,0,0,0) 
there. Let us then suppose that in the rest frame the 
flow velocity u^{x) is transformed in the following way: 



\x) h-> v ft (x ,x) = (j v ,j v v). 



(33) 



We use here the Minkowski metric instead of the general 
one considered so far because it can be approximated 
with the flat Minkowski metric ggf = diag(l, —1,-1,-1) 
for space-time sufficiently close to the space-time point 
x. If the LRC introduced into perfect fluid is isotropic, 
the v is also isotropic and v(x) = v(\x\), in the rest frame 
of U^(x). Notice that v(x) fully characterizes the LRC 
used. In the case considered here the local flow vector 



field u^{x) can be written as 

U^{x) = [7„,U|| +u ± ] 
with 7„ = u° = yj\ — + u±\ 2 and 



v\\+U 
IX 1 1 = — — — , u± 



1 + U ■ v\\ 



7c/ [1 + U ■ v\\] 



(34) 



(35) 



Here, v\\ = v\\ and v± = v± e± are, respectively, the 
parallel and perpendicular components of v with respect 
to U whereas e|| , ej_ are their unit vectors: 



"II = 7 =f ' V ^ = \A\I 1 ~ [iii • ey] 2 , (36a) 



and 



lu 



U, 



e i = 



(36b) 



Notice that the velocity field u^{x) (and, accordingly, 
also Su >1 (x)) are determined by the reference velocity field 
U ,J -(x) once we decide, by choosing v(x), on the type of 
the long-range correlations. As a possible example of the 
field v(x), let us consider 



v(x) = I(x) exp 



L(xY 



(37) 



where I{x) and L{x) are, respectively, the local correla- 
tion intensity and correlation length. In what follows we 
shall assume that only I(x) varies with the expanding 
fluid whereas the L is maintained as a global constant 
over the whole fluid evolution. In other words, we assume 
the intensity I(x) is determined by the dissipative fluid 
dynamics whereas the correlation length L is given by 
the particles composing the fluid. This assumption fits 
nicely the ridge phenomenon found at RHIC discussed 
in Section I, which indicates the persistency of the 
long-range correlations during the fluid evolution. 

To summarize this part, in our model its main quan- 
tity, Su^, is determined by only one equation Eq. (|30a[) 
(or equivalently Eq. (|3"Tj) ) and can be fixed by just one 
variable, the correlation intensity I(x). It means there- 
fore that correlation intensity itself is controlled by these 
equations and that it could be identified with a kind of 
"relaxation" to the equilibrium state taking place in the 
presence of LRC. 



C. Charge conservation 

Once <5u M is fixed, n and A are given by Eq. tfTO]) and 
(f!?3")) . respectively. As for Sn, it is determined by the 
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conservation of the (baryonic) charge. For the conserved 
charge current iV M , we can write explicitly that 



[{n^ + SnW + Sft")]. = 0, 



(38a) 
(38b) 



From the above two equations it is possible to write 



SN* = 0, 



(39) 



where SN 11 = N 11 — N£ q . Hence, the correction of the 
(baryonic) charge density in the presence of LRC is 
given by Eq. ([3T))) . 



V. OFF-EQUILIBRIUM ENTROPY CURRENT 

For the equilibrium thermodynamical quantities one 
has following fundamental thermodynamic relation, 

Ts cq = e cq + P cq - fJ-n cq , (40) 

which, using Eq. {JTJ) and (H~9|) . one can be rewritten as 



s eq (T,fJ,) = 



e + P — fin 
T 



411 - uSn 



T 



(41) 



Notice that first term at the r.h.s. of Eq. (|4Tj) can be re- 
garded as a natural extension of the equilibrium entropy 
density. We assume therefore that 



and, accordingly, that 



Ss 



e + P — [in 



411 - nSn 



(42a) 



(42b) 



Although there is no principle to determine the explicit 
form of the entropy flux $ M , a possible natural and simple 
candidate is 



Using it the explicit form of the off equilibrium entropy 
current is given by 

S» = {s cq + 6s)u» + 

= (1 - 7)[se q ^1 - «7[«c q t/ M ] - a[8N»] + U^, (44) 

where a = u/T and = j3u^ with (3 = 1/T. It 
should be noted here that the first three terms of Eq. (|4"4")l 
are quantities of the first order in 7, whereas the last 
term contains terms j 2 (see Eq. (|2"3")l and recall that 
-) = —U^Su* 1 is the first order infinitesimal quantity cor- 
rection induced by the LRC). It means that our model, 
in spite of applying simple form of the entropy flux "J^, 
introduces in a natural way the 2 nd order terms, both in 



the energy-momentum tensor, T^ v , and in the entropy 
current, 5 M . It is interesting to notice that 7 can be 
related to a derivative of the fluid velocity field u M , 



5u» ~ u£A", 



(45) 



where \^(x) is a kind of the local space-time scale. 
It connects Su^, which is crucial and most important 
parameter in our model, with the derivatives of the flow 
vector u^. 



The second law of thermodynamics, Sf^ > 0, can be 
expressed in our model as; 

[n/n-/i > + (d„a)[ 7 N^ + 6N»], (46) 

where current conservations, i.e., the second and third 
equations of Eqs.© and Eq. (f3"9")) are used. In the baryon 
free limit, one can also express the second law of ther- 
modynamics as an inequality of two kinds of time-like 
derivatives. Denoting ^ = u^X.^ and 9 = ut 1 ^, one may 
write 



n 



> 



2 1 



n 



n 



(47) 



where l/r^ = 9 - yfp The l.h.s. of Eq.(g7|| contains 
comoving time-like derivative d/dr with respect to the 
flow whereas the r.h.s. contains comoving time-like 
derivative with respect to the flow (denoted by dot). 
Note the minus sign in the Eq. ([4"T]) which implies that in 
the case of weak LRC (i.e., when I7I < 1), the inequality 
(second law of the thermodynamics) holds providing that 



n , dn n 

n + — > and — + — > 0. 



dr 



(48) 



This requirement, which is the consequence of our model 
presented here, restricts the relaxation process, namely 
the bulk pressure II should decrease with proper time r 
slower than exp[— t/t^] in the rest frame of and slower 
than exp[— r/rn] in the rest flame of Z7 M . 



VI. SUMMARY AND CONCLUDING REMARK 

We have discussed dissipative effects appearing by 
introducing the long-range correlations (LRC) between 
particles composing a perfect fluid. They arise be- 
cause the perfectness of the fluid is broken in such 
case, therefore such fluid obeys in general equations of 
imperfect hydrodynamics. In particular, to clarify the 
effects caused by the LRC, we have considered a thought 
experiment in which LRC can be switch 'on' and 'off'. 
It is assumed that the complete set of solutions for 
the perfect fluid (called RFF) obtained when LRC are 
not present are known. Then, the differences from the 
RFF, in velocity field, (5u M , in pressure, bulk pressure II, 
and in energy density, A, can be regarded as resulting 
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entirely from dissipative effects caused by the LRC. 

By switching off the LRC, the imperfect fluid should 
relax into a perfect fluid. Such a relaxation process 
can be regarded as a transformation of the energy- 
momentum tensor. It is assumed to be given by Eq.(|18p 
and we have shown that all dissipative terms in Eq.([2]), 
such as II, IT M and 7r Miy , are scaled by Su^. Moreover, 
one of dissipative equations Eq. ((la)) can be expressed by 
Ea.([5T]) (or, equivalently, by Eq . (|50a|l ) . 

It is noteworthy to observe that tensor relations 
Eqs. (|27f obtained in the presence of LRC, are the same 
as those found in the dissipative fluid corresponding 
to the non-extensive perfect fluid (g-fluid) discussed in 
Ref.gJ. 

When the LRC are spatially isotropic (in the rest frame 
of U^) with a constant correlation length L, we have 
shown that the local correlation intensity I(x) can be de- 
termined by Eq. (|31[) which is equivalent to dissipative hy- 
drodynamical equation Eq. (|la[) . Interestingly, we found 
that it is given by the form of constitutive-like equation. 
Moreover, it can also be presented in our model as expres- 
sion for the second law of the thermodynamics, Eq. (|47l) . 



The check of the causality and stability of solutions ob- 
tained in our model as well as numerical solutions are 
beyond the scope of this brief paper and will be discussed 
elsewhere. 

We close with observation that our model can be 
applied to a fluid with LRC originated from gravitational 
interaction. The perfect hydrodynamics applied for the 
universe in the early stage should be, for example, 
modified by such effects. However, this problem should 
be considered using General Relativistic (GR) hydrody- 
namics, as proposed in [5lM53l |. The GR hydrodynamics 
includes long-range interactions caused by gravity in 
a natural way. If LRC originated by the gravity also 
cause dissipative effects in a perfect fluid, then this 
may drive us to the interesting question whether GR 
hydrodynamics have certain correspondence in the 
(special) relativistic dissipative hydrodynamics. 
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